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Abstract 

We consider optimal stopping problems for a Brownian motion 
and a geometric Brownian motion with a "disorder", assuming that 
the moment of a disorder is uniformly distributed on a finite interval, 
f— | Optimal stopping rules are found as the first hitting times of some 

C/3 Markov process (the Shiryaev-Roberts statistic) to time-dependent 

<—| boundaries, which are characterized by certain Volterra integral equa- 

tions. 

The problems considered are related to mathematical finance and 
!— 1 can be applied in questions of choosing the optimal time to sell an 

I asset with changing trend. 

> 

Q\ Keywords: optimal stopping problems; disorder detection prob- 

lems; Shiryaev-Roberts statistic; Volterra equations. 

cn 

1 Introduction 

1.1. In this paper we consider problems of stopping a Brownian motion or 
a geometric Brownian motion optimally on a finite time interval when it has 
^ a disorder, i. e. its drift coefficient changes at some unknown moment of time 

from a positive value to a negative one. We look for the stopping time that 
maximizes the expected value of the stopped process. 

Let B = (B t ) t > be a standard Brownian motion denned on a probability 
space (Q, JP, P). Suppose we sequentially observe the process X = (X t ) t >o, 

X t = ntf + (n 2 - /ii)(t - 6)+ + aB t , 
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or, equivalently in stochastic differentials, 

dX t = [^I(t <9) + /i 2 I(t > 9)] dt + odB u X = 0, 

where p\ > > fi2, & > are known numbers and 9 is an unknown time of 
a disorder - a moment when the drift coefficient of X changes from value 
Pi to value fi2- The process X is called a (linear) Brownian motion with a 
disorder. 

Adopting the Bayesian approach, we let 9 be a random variable defined 
on (Q, JP, P) and independent of B. In this paper, in view of applications (see 
below), we assume that 9 is uniformly distributed on a finite interval [0,T], 
possibly, with mass at t — 0, T, i. e. the distribution function G(t) = P{9 < t) 
is given by 

G(t) = G(0) + pt for < t < T, G(T) = 1, 

where G(0) G [0, 1) is the probability that the disorder presents from the 
beginning, and < p < (1 — G(0))/T is the density of 9. The probability 
P(6> = 1) = 1 — G(T—) may be strictly positive. 

Let Tlx denote the class of all stopping times r < T of the process X. We 
consider the following two optimal stopping problems for X and the exponent 
of X (a geometric Brownian motion with a disorder): 

V®= sup EX r , V i9) = sup Eexp(X T -a 2 r/2). (1) 

The problems consist in finding the values V^ l \ and finding the stopping 
times t* \ ri^ at which the suprema are attained (if such stopping times 
exist). The superscript (/) stands for the problem for a linear Brownian 
motion, while (g) stands for a geometric Brownian motion. 

Observe that, roughly speaking, the processes X t and exp(X t — a 2 t/2) 
increase on average "up to time 9" and decrease on average "after time 9" . 
But since 9 is not a stopping time, we cannot simply take = 9 and need 
to stop by detecting the disorder based on sequential observation of X. 

We provide solutions to problems ([T]) using the results obtained in the 
recent paper (TO]. The central idea is based on a reduction to a Markovian 
optimal stopping problem using a change of measure. This approach has 
already been applied in the literature, but we were able to generalize it to 
the case of a finite time interval. 

1.2. For economic applications, the problems considered are related to the 
question when to quit a financial "bubble". By a bubble we mean growth 
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of an asset price based mainly on the expectation of higher future price; 
eventually a bubble bursts and the price starts to decline. 

Suppose that an asset price is modelled by a geometric Brownian motion 
with a disorder S = (S t )t>o- 

S t = exp(X t -aH/2), 

or, equivalently, 

dS t = ln x l{t <9) + fi 2 l(t > 9)] S t dt + aS t dB t , S = 1. 

Thus the price initially has a positive trend, which changes to a negative one 
at an unknown time 9. 

Let t = correspond to the "current" moment of time, when one is in 
a long position on an asset with positive trend. Usually, it is possible to 
predict that the trend will become negative by some (maybe, distant) time 
T in the future. Then one is interested in the question when it is optimal to 
sell the asset maximizing the gain. 

If nothing is known about the actual distribution of 9, it is natural to as- 
sume that 9 is uniformly distributed on [0,T] (since the uniform distribution 
has the maximum entropy on a finite interval). Interpreting the quantity 
ES T as the average gain achieved by selling the asset at time r, problem Q 
for a geometric Brownian motion seeks for the optimal time to sell the asset. 
The problem for a linear Brownian motion can be thought of as a problem 
of finding the optimal time to sell the asset provided that a trader has the 
logarithmic utility function, i.e. maximizes Elog(SV), which is equivalent to 
maximizing EX T with \j! x — H\ — u 2 /2, // 2 = // 2 — cr 2 /2. 

An interesting result that follows from the solution of problems ([TJ is 
the qualitative difference between risk-neutral traders who maximize ES T 
and risk-averse traders who maximize Elog(SV): if P(9 < 1) = 1 (i.e. the 
distribution of 9 has no mass at t — T), then a risk-neutral trader will sell 
the asset strictly before time T with probability one (P(ri 9 ' ) < T) = 1), while 
a risk-averse trader will wait until the end of the time interval with positive 
probability (P(rJ° = T) > 0); see the Theorem and Remark 1 in Section 

1.3. Problems Q was considered in the papers [U El E], assuming that 
the moment of disorder 9 is exponentially distributed. In PQ, the problem 
for a geometric Brownian motion was solved. It was shown that if /ii,/i2,cx 
satisfy some relation, the optimal stopping time is the first hitting time of the 
posterior probability process n t = P{9 < t \ J^ x ), where = a(X s ; s < t), 
to some level. In [2] this result was extended to all values of // 2 , c and the 
optimal stopping level was found. In [7j, the problem for a linear Brownian 
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motion was considered on a finite interval, i.e. assuming that one should 
choose r not exceeding some time horizon T (but the disorder may happen 
after T). It turned out that the problem is equivalent to the original Bayesian 
setting of the disorder detection problem when one seeks for a stopping time 
minimizing the average detection delay and the probability of a false alarm 
(see e. g. [H IH [9]). The paper [7] also briefly discusses the optimal stopping 
problem for a geometric Brownian motion on a finite interval, but does not 
provide an explicit solution. 



2 The main result 

2.1. Let fi = (fii — [i^jo denote the signal-to-noise ratio. For convenience 
of notation, introduce the process X = (X t ) t > , X t = (X t — nit)/a, which is 
a Brownian motion with the unit diffusion coefficient and the drift coefficient 
changing at time 9 from value to value (— /x). 

Introduce the Shiryaev-Roberts statistit^ip = (ipt)t>o' 

ip t = exp(-/jX t - /iH/2) (^o + pJ q ex P {yX s + ^ V 2 ) ds^j (2) 

with ip = G(0). Applying the Ito formula it is easy to see that ip satisfies 
the stochastic differential equation 

dip t = pdt — fj,ip t dX t . (3) 

On the measurable space (Q, J^jO, = a (-^-t] t < T), define the prob- 
ability measures PW and P (s) such that X t is a standard Brownian motion 
under P® and (X t — at) is a standard Brownian motion under P^. These 
measures will be used to solve the problems and respectively. It 
is well-known (see, e.g., Ch. 7]) that PW and are equivalent on the 
space (fi, ^t)- 

For any x > 0, by Ei [•] and Ei [•] we denote the mathematical expec- 
tations of functionals of the process (iftt)t>o defined by (j2])-([3]) with the initial 
condition ipo = x, when X is respectively a standard Brownian motion or a 
Brownian motion with drift a. For brevity, instead of EgLj-] and E^L[-] 

we simply write E^ [ ■ ] and E^ [■]. 

The main result of the paper is the following theorem. 



In a general case, the Shiryaev-Roberts statistic is given by ipt ~ 
{dP^/dPf)Jg(dPf/dP° s )dG(s), where P| = P s | J^ x for the measures P° or P°°, 
corresponding to that the disorder happens at time t = or does not happen at all (for 
details, see [10] and the proof of the Lemma in Section [3j). 
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Theorem. The optimal stopping times in the problems V® and V^ 9 ' are 
given respectively by 



r. 



T 



W = inf{t > : if) t > a (,) (t)}AT, 
(9) = inf{t > : ip t > a {9 \t)} A T, 



where a®(t), a {9 \t) are non-increasing functions on [0, T] being the unique 
solutions of the integral equations (t e [0, T]) 



J E$ )(t) [(//! - (/X! - /i 2 )^)l{^ s < a®(t + = 0, (4) 

/ T i ES )(t) [e^( yUl (l-G(t + S ))-| yU2 |^)l{^<a^(t + S )}]rf S = (5) 
in £/ie c/ass of continuous bounded functions on [0, T) satisfying the conditions 
a {l) (t) > forte [0,T), a«(T) = , (6) 

a^(t)>^(l-G(t)) forte [0,T), a W CO = /^(l - G(T-)). (7) 

| A*2 | | A*2 | 

The values V® and V® can be found by the formulas 

V® = f T ^[( m - - fM 2 )if) s )l{if) s < aV\s)}]ds, (8) 
Jo 

V® = 1+1 E^[e^(Ml - G{s)) - |/i 2 |^)IR < a<*)(s)}]ds. (9) 
Jo 

Remarks. 1. Regarding the difference between a risk-averse trader and 
a risk-neutral trader mentioned in Section [TJ observe that if G(T—) = 1, 
then P(ri 9) < T) = 1 because a^(T) = 0, while P(rJ° < T) < 1 because 
a^(T) > (in the latter case the process ip stays below a^(t) on the whole 
interval [0, T] with positive probability). 

2. In the above mentioned papers [UEIII], solutions to problems ([!]) when 
9 is exponentially distributed were given in terms of the posterior probability 
process ir t = P(9 < t | Using the Bayes formula, one can check that 

the processes 7r and if) are connected by the formula if> t = 7it(l — G(t))/(1 — Tr t ) 
(see [10J). Consequently, it is easy to reformulate the Theorem in a such way 
that tP and are the first moments of time when the process re crosses 
time-dependent levels. We prefer to work with the process if) because it has 
a somewhat simpler form than it, when 9 is uniformly distributed. 
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2.2. Equations (J4J), ([5]) can be solved numerically by backward induction: we 
fix a partition < to < ti < • • • < t n = T of [0, T] and sequentially find the 
values a(t n ), a(t„_i), . . . , a(t ). The value a(T) can be found from condition 
^ or ([?]) respectively. Having found the values a(tfc), a(£fc+i), • • • , a(t n ) and 
numerically computing integral Q or ([5]) for t = tk-i through the values of 
the integrand at points tk-i,tk, ■ ■ ■ ,t n , we obtain the equation, from which 
the value a(t fc _ 1 ) can be found. Repeating this procedure, we find the value 
of a(t) at every point of the partition. 

To compute the mathematical expectations in Q, ([5]), ([8]), rt9j), one can 
use the Monte-Carlo method or use the explicit formula for the transitional 
density of ip (see e.g. [10J, where it was obtained from the joint law of an 
exponent of a Brownian motion and its integral that can be found in [6]). 

As an example of a numerical solution of equations Q, (pi), on Figure [l] 
we present the optimal stopping boundaries for the case T — 1, //1 = 1, 
^2 = -1 ; a = 1, G(t) = t for t < 1. 




3 Proof of the theorem 

3.1. To prove the theorem, we first reduce problems and to optimal 
stopping problems for the process if>, and then apply the Proposition from 
the Appendix, which was proved in [10] . The method we use is based on the 
ideas of [U E] . 
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Lemma. The following formulae hold: 

V (l) = sup E {1) f (//i - (/ii - fi 2 )ip s )ds, 
tsOTt Jo 

V {9) = 1 + sup / e w '[/ii(l - G(s)) - \fi2\ips]ds. 
tgW-t Jo 

The supremum in each formula is attained at the stopping time which is 
optimal in the corresponding problem ([I]). 

Proof. It is sufficient to show that for any stopping time r G OTTy 

EX T = E« I (/n - (m - ^)ijj s )ds, (10) 
Jo 

Eexp(X T -a 2 r/2) = 1 + E (s) / e^ s [^(1 - G(s)) - |/x 2 |Vj ds. (11) 

On the measurable space (fl, ) define the family of probability mea- 
sures (P")o<m<t such that under P" the disorder occurs at the fixed time 
u, i. e. for each < u < T the process X can be represented as X t = 
Hit + (/i2 — Hi)(t — u) + + cB t , where is a standard Brownian motion 
under P". By E u [-} we denote the mathematical expectation with respect to 
P" and by P 4 , P", Pf\ P[ g ^ we denote the restrictions of the corresponding 
measures to the a- algebra = cr(X s ; s<t),0<t<T. 

Let us prove (10). Since the Brownian motion B is a zero- mean martin- 
gale, we have 

EX T = E [^r + (/i 2 - /i!)(r - 9) + ] . (12) 
Consider the second term in the sum: 

E(r-0) + = / E u [{T-u)l{r>u)\dG{u). (13) 
Jo 

Observe that for any < u < T the following relation is valid: 
E u [{t-u)I{t >u)\ = j E u I(s<r)ds 

J u 

= f E w [ex P (-n(X s -X u )- fi 2 (s-u)/2)l(s<r)]ds, (14) 

J u 

where we use that I(s < r) is an ^^-measurable random variable and 

dP u ~ ~ 

- exp(-^(X s -X u )-/i 2 ( S -n)/2) 



dP ( P 
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(the explicit formula for the density of the measure generated by one Ito 
process with respect to the measure generated by another Ito process can be 
found in e. g. [3]). From (13)-(14), changing the order of integration we find 



T - 



E (0 



JO 



exp(— fi(X s — X u ) — ji 2 (s — u)/2)dG(u)ds 



E® / fads. 



Then using (12), we obtain representation (10). 
Let us prove (11). We have 



dP't 
dP[ a) 



exp(-aX t + a 2 t/2) ■ exp(-fi(X t - X s ) - fi 2 (t - s)/2), s < t, 
exp(-aX t + a 2 t/2), s > t, 



which implies 

dP t r T dp: 



dp 



(.</) 



o dP\ 9) 



dG{s) 



exp ( - aX t + a 2 t/ 2) 

exp(-/i(X t - X s ) - ^{t - s)/2) dG{s) + 1 - G(t) 
= exp(-aX t + a 2 t/2) + 1 - G(t)). 
Consequently, 

Eexp(X T - a 2 r/2) = E« exp(X T - a 2 r/2)\ 

= E^ [e» 1T (^ T + l -G(t))]. 
Applying the Ito formula, we get 

e" 1T (f + 1 - G(t)) 

= 1+1 e ws [// 2 -0 s + - G(s))] ds + ji I if; s d(X s - as). 



Taking the mathematical expectation E^[-] of the both sides and using 
that (X s — as) is a standard Brownian motion under P*^ and, hence, the 
expectation of the stochastic integral equals zero, we obtain (11). □ 
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3.2. Now the proof of the Theorem follows from the Proposition in the 
Appendix - for the problem V® we use that the process ip satisfies equation 
@ with X being a Brownian motion under the measure P"', and for the 
problem we use that ip satisfies the equation 

dip t = (p- (/ii - p 2 )'4 , t)dt - pip t d{X t - at), (15) 
where (X t — at) is a Brownian motion under P^ 9 \ 
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Appendix 



Let B = (B t ) t >o be a standard Brownian motion defined on a filtered prob- 
ability space (fi, J£~, (J^t)t>o, P) and ip = (ipt)t>o be a stochastic process sat- 
isfying the stochastic differential equation (cf. (15)) 



dipt — (p + bip t )dt — pip t dB t 



(16) 



where b, p, G M and p > 0. 

Consider the optimal stopping problem consisting in finding the quantity 



V 



sup E 

T<T 



(f(s) ~ ^s)ds, 



(17) 



where A G K, and f(s) is a non-increasing bounded function which is contin- 



uous and strictly positive on [0,T). The supremum in (17) is taken over all 
stopping times r of the filtration (<^t)t>o satisfying the condition r < T a.s. 

For arbitrary x > 0, let E x [-] denote the mathematical expectation of 
functionals of the process (ipt)t>o, satisfying (16) with ip Q 

The following result was proved in |10j . 



x. 



Proposition. The optimal stopping time in problem (17) is given by 

r* = inf{t >0:ip t > a(t)} A T, 
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where ait) is a non-increasing function on [0,T] being the unique solution of 
the equation (t G [0,T]^) 

T-t 

E o(t ) [e Xs (f(t + s)- ^ s )l{^s < a(t + s)}] ds = 

m the class of continuous bounded functions on [0, T) satisfying the conditions 

a(t)>f(t) forte [0,T), a(T) = /(T-). 
T/ie quantity V can be found by the formula 

V= f E A) [e Xs (^ s -f(s))I{^ s <a(s)}}ds. 
Jo 
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